Finite-size corrections of an integrable chain with 
alternating spins 



B - D D5rfelt and St MeiBner§ 

Institut fiir Physik, Humboldt-Universitat , Theorie der Elementarteilchen 
Invalidenstrafie 110, 10115 Berlin, Germany 

Abstract. In this paper we calculate the finite-size corrections of an anisotropic 
integrable spin chain, consisting of spins s = 1 and s = ^ . The calculations are done 
in two regions of the phase diagram with respect to the two couplings c and c. In case 
of conformal invariance we obtain the final answer for the ground state and its lowest 
excitations, which generalizes earlier results. 



PACS numbers: 75.10 JM, 75.40 Fa 



Short title: Finite-size corrections 
February 7, 2008 



X E-mail: doerfel@qft2.physik.hu-berlin.de 
§ E-mail: meissner@qft2.physik.hu-berlin.de 



2 

1. Introduction 



In 1992 de Vega and Woynarovich constructed the first example of a spin chain with 
alternating spins of the values s = ^ and s = 1 on the basis of the well-known XXZ(|) 
model. We call this model XXZ(i, 1). Later on a lot of interesting generalizations have 
been presented p, |, ^. After de Vega et al p, |] we have studied the XXZ(|, 1) model 
in two subsequent publications 0, §]. In our last paper p we determined the ground 
state for different values of the two couplings c and c (for the details see section 3 of that 
paper). Disregarding two singular lines we have found four regions in the (c, c)-plane 
which can be divided into two classes. The division is made with respect to the occurance 
of finite Fermi zones for Bethe ansatz roots. The two regions with infinite fermi zones 
only are well studied in the framework of Bethe ansatz. On that basis we consider 
the finite-size corrections for the ground state and its lowest excitations using standard 
techniques 0, |ll|. It is remarkable that they allow to obtain an explicit answer only 
in the conformally invariant cases, which are contained in the two regions considered. 
The results can easily be compared with the predictions of conformal invariance. 

The paper is organized as follows. Definitions are reviewed in section 2. In section 
3 we calculate the finite-size corrections for both couplings negative. The same is done 
in section 4 for positive couplings. Here it was necessary to set c = c to obtain explicit 
answers. Section 5 contains the interpretation of the results and our conclusions. 



2. Description of the model 

We consider the Hamiltonian of a spin chain of length 2A^ with even 

n{-f) = cH{-f) + cn{-f). (2.1) 

The two Hamiltonians can (implicitly) be found in paper [|l|, they both contain a two- 
site and a three-site coupling part. Their explicit expressions are rather lengthy and do 
not provide any further insights. They include a XXZ-type anisotropy parametrized by 
e*''', we restrict ourselves to < 7 < 7r/2. The isotropic limit XXX(|, 1) is contained 
in 0. The two real coupling constants c and c dominate the qualitative behaviour of 
the model. The interaction favours antiparallel orientation of spins, for equal signs of 
the couplings its character resembles ordinary XX Z model. A new kind of competition 
comes in for different signs of couplings where the ground state is still singlet but with 
a much more involved structure. 

The Bethe ansatz equations (BAE) determining the solution of the model are 

/ sinh(Aj + q) sinh(Aj + 27) \ ^ ^ -rj sinh(Aj - + 27) ^ 
\^sinh(Aj — ij) sinh(Aj — ^7) J ^J--^ sinh(Aj — A^ — 27) ' 
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One can express energy, momentum and spin projection in terms of BAE roots A^: 
E = cE + dE, 

2 sin 7 



' 5 

M 



E = -y- 

cosh 2Xj — cos 7 

. (2.3) 

jT^ cosh 2Xj — cos 27 

P='-f L + '"> U lu ( "°"<'^ + '^> H (2 4) 

2,ta"'Umh(A,-.i)j+"'Uw'(Aj-n)jr * ' 

5. = ^ - M. (2.5) 

We have defined energy and momentum to vanish for the ferromagnetic state. The 
momentum operator was chosen to be half of the logarithm of the 2-site shift operator 
which is consistent with taking the length of the system as 2N instead of A^. 

3. Calculation of finite-size corrections for negative couplings 

In section 3 of paper |Q] we have carried out a detailed analysis of the thermodynamic 
Bethe ansatz equations (TBAE) at zero temperature and obtained the ground state. 

We found a large antiferromagnetic region in the (c, c)-plane (depending on 7) where 
the ground state is formed by roots with imaginary parts |, the so-called (1, — ) strings. 
The Fourier transform of their density is given by [0 

l + 2cosh(p7/2) 

Po{P) = 7; r^T w7^- 3.1 

2 cosh(p(7r — 7)/2) 

Depending on the signs of c and c the region is described by the connection of three 
parts: 



a) 



b) 



£ < 0,c < 0; 



c < 0, c> 



c 1 , 2tt 

— < for < 7 < — 

\c\ 2 cos 7 5 
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c) 

c > 0, c < 



IcI 8cos^7 ^ TT 

V > ^ for < 7 < 77 

c 4 cos^ 7 — 1 3 

|c| 2 ^ TT TT , , 

U> ^ for 77<7<7T- 3.2 

c COS 7 3 / 

Here for shortness we have introduced 

7 = ./^ V (3.3) 
2(7r -7) 

We shall now calculate the finite-size corrections for the ground state and its excitations. 
In paper the structure of excitations in the framework of BAE roots was obtained 
for c < 0, c < 0. Our results immediately apply to the whole region (3.2), because we 
had to ensure only that the ground state consists of (1, — ) strings which follows from 
TBAE. 

Because we are interested in the lowest excitations only, we disregard the bound 
states 0] and consider those excitations given by holes in the ground-state distribution 
which are located right (or left) from the real parts of all roots. The number of those 
holes we call {H~). We follow the standard techniques developed in and |[T0[| . 

For transparency we employ the notations of lll|] as much as possible. We decompose 



(^N = p'i^ + p'i^ + A(T,v (3.4) 

where the upper index describes the two terms on the RHS of formula (|3.1| ). The basic 
equations are then 

^^"^-6^ = 077 |_^t/Apf,^)(A)|l^5(A-Afc)-a^(A)| 



2N 



-ax 



d\pf^ W[]^Y. ^(^ - ^k) - on{X) I (3.5) 



and 



(3.6) 



Aa^(A) = - d^^p{X - p) |l ^ 5(A - A,) - ajv(A)| 

where the Fourier transform of the kernel p{X) is given by 

p/ N ^ sinh((7r/2-7)^) 

^ ' 2sinh(cJ7/2)cosh(cj(7r/2-7/2))" ^ ' 

The summation on the RHS of (3.5) and ( p. 61 ) is carried out over the real parts of all 
the roots (without the holes). Using Euler-MacLaurin formula equations (3.5) and ( |3.6| ) 
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are rewritten as usual, e.g. 
a^(A)-p«(A)-pf)(A) = 

1 1 



Loo ^f^^^^^^'fP^^ - ^) - ^P(^ ~ 



(T7v(A 



p{\ - A+) 
—p{X - A- 



(3.1 



Here (A^) is the real part of the largest (smallest) root. For A > A+ the part in 
round brackets can be omitted and equation (3.8) after a shift converts into a standard 
Wiener-Hopf problem to be solved. In the expression for AEn both parts have to be 
kept, so we need the solution for A > A"*" and A < A~, which are simply related (but 
not equal) by symmetry. 

For the solution with A > A^ we define as usual 



e^^V^(A + A+)rfA, 



(3.9) 



a^(A + A" 



ajv(A + A+) for A > 

< 
> 







for A 








(3.10) 



After Fourier transformation equation (3.8) takes the form 

X_{uj) + (1 - P(cu))(X+(cu) - Ciu)) = F+{u) + F^u) - C{u) (3.11) 

where we have marked all given functions of our problem by a bar. F±{uj) are defined 
as above using instead of aN the sum p'q^ + Pq'\ Further 

1 iuo 



^^'^^ 2N ' 12iW^(A+) 
Now we have to factorize the kernel 



(3.12) 



(3.13) 



with G±{uj) holomorphic and continuous in the upper and lower half-plane respectively. 
Noticing that 



P{uj,'j) = K{uj,n- 7) 



(3.14) 



where K is the analogous function in paper [|TT| we take the factorization from there. 



2^r ( 1 - e^(-) 



^ / 1 _ i(7r -7)^ \ ^ n i-fu 

I 2 277 / V2 277 



G_{-u), (3.15) 



Y 



In 



77 



.7, 



77 — 7 



In 



77 



77 — 7 
. 7 . 



(3.16) 
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It is chosen to fulfill 
where 

91 



~ i H 



+ — 



I 

12 



7 . 



vr — 7 

After the necessary decomposition 

G-{uj)F+{uj) = Q+{uj) + Q.{uj) 
equation (3.11) has the desired form 



Q^{lo) = Q_{lo)-G-{lo) X^{lo) + C{uj)-F^{uj) = P{lo) (3.20) 



leading to an entire function P{uj) given by its asymptotics. 

p/ X ^ m \ 

^'^^ 12A^Vjv(A+) 2N 12ArViv(A+)' 

Equation (3.20) yields the solution for X^{uj): 

X+{uo) = Ciu) + G+iu) [p{lj) + Q+{u) . 
for our purposes it is sufficient to put 



and hence 



0+(^) 



3 7rA+/(7r— 7) 



'1 + 2 cos 7) 



TT — iuj{n — 7) 
G+(z7r/(7r - 7))e-"^+/("-^)(l + 2 cos 7) 



TT — iLu{TT — 7) 

Next we must determine by normalization the value of the integral 

CTAr(A)(iA = Zn{oo) — Z]\f{A'^). 



A+ 



After a thorough analysis we found for our case that the relation 



2" ~ 



uS 



1 



(3.17) 



(3.18) 



(3.19) 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



holds where u = ^. Nevertheless we shall not claim formula ( p.25|) for all possible states 
T0| , especially we expect effects like those described in [|12| for higher excitations. Then 



we have 



2;^(±oo) - 2;;v(A^) = ± 



N \ 2 



1 c: . 



H 



V2 



(3.26) 



yielding the important equation 

G+(m/(7r - 7))e-'^^^/("-^) (1 + 2 cos 7) 



TT 2N 12A^V^(A+) N 

The other normahzation equation is obviously 



(3.27) 



(3.28) 



Now we can proceed in the usual way keeping in mind the changes arising especially 
from the last two equations. 



N 



TT 



2N 



TT — 7 



(c + 2ccos7)G'^ 



ITT 



TT — 7 



P 



TT — 7 



ITT 



TT — 7 



-7rA+/(7r— 7) 



After some algebra using equations ( p.27| ) and (|3.28|) this turns out as 
TT^ (c + 2ccos7) 



(3.29) 



AE 



N 



2N 



TC 



X 



7 1 + 2 cos 7 
1 

+ 



X 



+ 



+ 



1 



Finally, for further interpretation we put it in the form 

A^AT 2c + 4c cos 7 



(A- 



TT 



2A^ 



1 27r 

l + 2cos7 7r-7l Q AN^ ^ AfP 



A2^ 
2 ^ 2z/ , 



(3.30) 



(3.31) 



with A = {H^ — H~')/2 as an integer number. 

The momentum correction AP/v is obtained from relation (3.5) after substituting 
the hole energy Sh = —2ctipq^^ — 2cnp^^^ by the hole momentum p/i(A) = 
I arctan(sinh(7rA)/(7r — 7)) + arctan(sinh(7rA)/(7r — '^))/ cos 7)+ const (see 0). 

Comparing the asymptotics for large A of both £/i(A) and p/i(A) gives the speed of 
sound and helps to shorten the calculation of AP^r. Therefore 



AP 



N 



2N 



TT J 1 

2 liiVV 



(A-)2 - (A^ 



const. 



(3.32) 



We are not interested in the constant term, being some multiple of tt. 
Finally 

27r 



AP 



N 



2N 



:S.A. 



(3.33) 



The interpretation of our result will be given in section 5. We stress once more, that 
to obtain equations (3.31) and ( |3.33| ) it was not necessary to put c = c. The coupling 
constants are only constrained to stay in region (3.2). 



4. Calculation of finite-size corrections for positive couplings 

Now we consider region c > 0, c > and rely on the analysis of paper 

The ground state is given by two densities, a^^^'^\X) for the real roots and cr^'*(A) 
for the real parts of the (2, +) strings. One has 

^^/^HA) = ag\X) = ^ = .(A). (4.1) 

27cosh(7rA/7j 

The physical excitations are holes in those distributions. As in section 3 we consider 
only holes situated right (or left) from all roots. With the usual technique and the 
results of |I| we have obtained after some lengthy but straightforward calculations the 
basic system for the density corrections 



/•oo ( 1 

^ E '^l-" - A.) - 
- df^ r(A - /x) 1 1 f: 6{f, - Q - a«(/i) | . (4.2) 



We have denoted the real roots by Aj (their number is M1/2) and the real parts of the 
strings by C,j (number Mi). The function r(A) is given via its Fourier transform 

R( ) = sinh(a;(7r - 3j)/2) 

2sinh(cj(7r-27)/2)cosh(cu7/2)" ^ ' 



The energy correction takes the form 

^ = - iyHx)U'f:H^ - A.) - .'^W 

-^i r d\s{\)l-^TS(\-i,)-a'^'{X)\. (4.4) 



Once more we shall follow |]TT[. The maximal (minimal) real roots we call Ay^ and for 
the strings we use Af respectively. Instead of one C{iu) we have now Ci{iu) and Ci/2(i^) 
generalized in an obvious way. The same applies to F{uj). The main mathematical 
problem is the factorization of a matrix kernel 

{1-K{uj)y^ = G+{u)G4uj) with 

G.{uj) =G+{-ujf (4.5) 



see 



I) and 







S{uj)e 



G+ is now a matrix function and G\ stands for its transposition. The two component 
vector Q^{oj) is (see equation ( |3.19|) ) 

7rA+/7 



(4.6) 



G+(z7r/7) 



vr — zci;7 \ e 

As usual we define the constant matrices G\ and G2 by 

1 



(4.7) 



and as before one has 



G2 — 77 G*!- 



1 + Gi- + G2-^ 



The two component vector P{uj) is then 
Piuj) 



/ 1 joj 

2N 



Gi 



(4.9) 



(4.10) 



12JV2<xW(A+) 



(4.11) 



12JV2<7(,i)(A+) 

and therefore the shifted densities are expressed in the form 

Now it is necessary to find the analogue of equation (3.26) for the two counting functions. 
Here it would be necessary to consider different cases depending on the fractions of 
uSz/N or 2vSz/N. From our experience we know that the result of the finite-size 
corrections does not depend on those fractions, while relations hke (|3.25|) obviously do. 
Being interested only in the former we shall proceed as easy as possible and consider 
only the case with vanishing fractions. 



,(1/2), 
-N ' 



±ool 



±^ - + ^V2 



(4.12) 



z«(±oo) - 4^(Af ) = ±^ - + H( 

Easy counting leads to expressions for the numbers of the holes 
Hi = 2Sz 

Hi/2 = 2Sz + 2Mi - N. 



(4.13) 



(4.14) 
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We expect their modifications for non-vanisliing fractions. We stress tliat botli numbers 
are even. 

Equation ( p.27| ) is now more complicated 



g-7rA+/7 



G-\0)B^+{ 2f 



with the definitions 



and 



1_ / -uS, + H^/, 



N 



-2vS, + Hf 



Ht — Hi 



2N 



1 

N 



12iVV(^>(A+) 



S,-iyS, + Mi-^± A(V2) 



(4.15) 



(4.16) 



(4.17) 



The other normalization condition is obviously 

1 

+ ^^1 



(1/2). A + 



N 



12) 




12JV2a(^/2)(A+ ) 



^ \ , G+(m/7)^ 



1 

2Af 



7 



-7rA+2/7 



e 1/2 
g-7rA+/7 



(4.18) 



After combining equations (4.18) and (4.15) we obtain from equation (4.4) 



27V ~ 7 I ge^^^/^ 




12Af2 



1 



\2N^a^^\kt) 



(A+ ^ A",fi+ ^ 5" 



(4.19) 



This result is valid for any positive c and c. As in paper |^ no further progress can be 
made unless the factorization is explicitly known which is not the case up to now. For 
c = c = c (conformal invariance) the problem simplifies and the final answer can be 
obtained. 

After some lengthy calculations one arrives at 



AE 



AT 



2A^ 



7 



12]V^ + H(^^^ 



5 



BtB 



[Btf 



TT — 37 

2(7r-27), 



[B' 



.B-) 
(4.20) 



which can be brought into the form 



AE 



N 



2N 



27TC 

7 



vr 2 27r 

+ 



6 4A^2 4JY2 



^{l-2u)Sl + ^{Hy,-^ 



4 



+ (A(V2))2 _ ^(1/2)^(1) ^ 1 (^{1))2_L_ 



(4.21) 
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For c 7^ c we expect the result to be much more comphcated but of the same order 

As above the momentum correction is given from relation (4.4) after replacing 
e^^^"^^ = 2cTTs{X) by p^^^"^^ = arctan e'^'**/'^ and e^^'* = 2c7rs(A) by pj^^'' = arctan e'^'*'/'*'. 
The values of the hole momenta are taken from where an additionel factor ^ must 
be introduced to take into account our definition of momenta (p^). 

As in section 3 comparing the asymptotics for c = c gives the speed of sound and 
together with equation (4.20) the momentum correction 
AP/v vr n 



2A^ 2 12 



'2 



+B^B^-BtBt-mf-{Bt^'^ " 



'2(7r - 27) 

Disregarding the constant (multiple of vr) we have 

Hi Hi/2 



+ const. (4.22) 



2 



(4.23) 



5. Conclusions 



In sections 3 and 4 we have determined the finite-size corrections of our model for 
two different cases. Equations (3.31) and (|3.33|) give the result for region (3.2), while 



equations (4.21) and (4.23) are valid for c = c = c> 0. In both cases we have conformal 
invariance. That seems to be the reason why during the calculations a lot of terms 
cancel each other and the result is considerably simple. 
From the asymptotics of Eh and ph we find 
2n C + 2CCOS7 

Vs = z- > (5.1) 

TT — 7 1 + 2 COS 7 

for the speed of sound and from relation (3.31) the value of 1 for the central charge. 
Formula ( ^.1|) generahzes the former result obtained for equal (but negative) couplings 
0. Analogously follows 

V. = ^-^ (5.2) 

7 

and the central charge equals 2 |]l|. 

For completeness we mention the heat capacities per site at low temperature given 

by 

C = ^ (5.3) 

where we have denoted the central charge by to avoid confusions. Therefore 

1 + 2COS7 (7r-7)T 
c + 2c cos 7 6 
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C = f (5.5) 



and 



in agreement with former results Formula (|5.4| ) generalizes our calculations for c = c. 

The dimensions Xn and the spins s„ of the primary operators follow from formulae 
(3.31) and ( |333|) . 



= :££i^ + (5.6) 

2 2z/ ^ ^ 

Sn = S,\A\ (5.7) 

for negative coupling. It is interesting to compare this with the result for the XXZ(|) 
model, where u is simply replaced hj 1 — u. Formulae ( |5.6|) and ( p.7|) have to be 
understood in the way that for general excited states (arbitrary holes and complex 
roots) Sz and A are replaced by more complicated integer numbers ||10|| . 
For positive couplings 

X. = \{1 - 2u)S! + \ - f )' + (A(^/^))' - A(V2) AW + 1 (A(^))^ (5.8) 

When z/ — > the first two terms agree with paper |^ while the other terms are connected 
with the asymmetry of the state which was not considered there. The dimension of 
a general primary operator depends on four integer numbers. The second of them 
measures an asymmetry between the number of holes among real roots or strings 
respectively. Once more, for more complicated states the integers in equation (5.8) are 
replaced by other ones depending on the concrete structure of the state. We mention 
that relation (5.8) can be " diagonalized" to resemble the expression of two models both 
of central charge 1. 

"^^"2 2 ^^ + 2 2 +2^ + 21-2^.^ ^^'^^ 

with 

2 ' 

A(i) 

A, = (5.10) 

Expression (5.9) becomes even more symmetric if one remembers the first equation of 
(4.14) and the definition ( [4.17| ). Then twice a certain number of holes is linked with its 
appropriate asymmetry. 

From formula (4.23) the spins of the primary operators are 



(5.11) 
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and after using (5.10) 

Sn= Ai ( Hi/2 



Hi 
2 



(5.12) 



with the same symmetry as relation (5.9). 

Finally, we shall determine the magnetic suszeptibilities per site at zero temperature 
and vanishing field from our finite-size results. That can be done, because the states we 
have considered include those with minimal energy for a given Sz (magnetization) [13 
Differentiating twice the energy with respect to Sz gives the inverse susceptibility. 

Hence 

TT — 7 / 1 + 2 COS 7 \ 11 
47r7 I c + 2c cos 7 / t's 27 



X 



(5.13) 



and 



1 



X = 7. TT- = ^ (5.14) 

2c7r TT — 27 t;<j TT — 27 

for the two cases respectively in agreement with earlier results 0, ^ H. Expression 
( ^.13| ) for general couplings in region (3.2) had not been derived before. 
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